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Abstract 



We study the phenomenology of a new Minimally-extended Supersymmetric Standard 
Model (nMSSM) where a gauge singlet superfield is added to the MSSM spectrum. 
The superpotential of this model contains no dimensionful parameters, thus solving 
the /x-problem of the MSSM. A global discrete i?-symmetry, forbidding the cubic 
singlet self-interaction, imposed on the complete theory, guarantees its stability with 
respect to generated higher-order tadpoles of the singlet and solves both the domain 
wall and Peccei-Quinn axion problems. We give the free parameters of the model 
and display some general constraints on them. A particular attention is devoted to 
the neutralino sector where a (quasi-pure) singlino appears to be always the LSP of 
the model, leading to additional cascades, involving the NLSP — > LSP transition, 
compared with the MSSM. We then present the upper bounds on the masses of the 
lightest and next-to-lightest - when the lightest is an invisible singlet - CP-even 
Higgs bosons, including the full one-loop and dominant two-loop corrections. These 
bounds are found to be much higher than the equivalent ones in the MSSM. Finally, 
we discuss some phenomenological implications for the Higgs sector of the nMSSM 
in Higgs production at future hadron colliders. 



1 Introduction 

Supersymmetry provides a well defined framework for the study of physics beyond the 
Standard Model (SM). Its main motivation has been the special properties of super- 
symmetric (susy) theories with respect to the hierarchy problem. In addition, the low 
energy data support unification of the gauge couplings in the susy case, in contrast to 
what happens in the SM scenario. Another interesting feature of susy models is that 
the breaking of the electroweak (EW) symmetry can be radiatively triggered by the 
largeness of the top quark mass [0] . The Minimal Supersymmetric Standard Model 
(MSSM) 0] is defined by promoting each standard field into a superfield, doubling 
the Higgs fields and imposing i?-parity conservation. Due to the non-observation of 
superpartners of the standard particles, supersymmetry has to be broken at a scale 
M susy not larger than O(TeV), so that it still provides a natural solution to the hi- 
erarchy problem. Unfortunately, a phenomenologically acceptable realization of EW 
symmetry breaking in the MSSM requires the presence of the so-called //-term, a 
direct susy mass term for the Higgs fields, with values of the (theoretically arbitrary) 
parameter /i close to M susy or Mw, when its natural value would be either or Mp. 
Of course, there exist explanations for an 0(My/) value of the //-term, alas, all in 
extended settings ||. 

The more or less straightforward solution to the //-problem is to promote the //- 
parameter into a field whose vacuum expectation value (v.e.v.) is determined, like 
the other scalar field v.e.v. 's, from the minimization of the scalar potential along the 
new direction || j|, |[ [7||. Naturally, it is expected to fall in the range of the other 
v.e.v. 's, i.e., of order 0(M susy ). Such a superfield has to be a singlet under the SM 
gauge group. In order to avoid introducing new scales into the model one should 
stick to dimensionless couplings at the renormalizable level. This can be achieved 
by imposing a ZZ 3 symmetry on the renormalizable part of the superpotential. The 
resulting model, the Next-to-Minimal Supersymmetric Standard Model (NMSSM), 
has the following superpotential: 

W = XSH 1 H 2 + ^S 3 + ... (1) 

where the dots stand for the usual quark and lepton Yukawa couplings (cf eq. (0)). 
The ZZ 3 symmetry is spontaneously broken at the EW scale when the Higgs fields get a 
non-zero v.e.v. It is well known, however, that the spontaneous breaking of such a dis- 
crete symmetry results in disastrous cosmological domain walls, unless this symmetry 
is explicitly broken by the non-renormalizable sector of the theory. Domain walls can 
be tolerated if there is a discrete-symmetry- violating contribution to the scalar poten- 
tial larger than the scale 0(1 MeV) set by nucleosynthesis ||. Heavy fields interacting 
with the standard light fields generate in the effective low-energy theory an infinite 
set of non-renormalizable operators of the light fields scaled by powers of the char- 
acteristic mass-scale of the heavy sector (Mp, Mqut,---)- These terms appear either 
as .D-terms in the Kahler potential or as F-terms in the superpotential. It is known 
however that gauge singlet superfields do not obey decoupling ||, [|, so that, when 
supersymmetry is either spontaneously or softly broken, in addition to the suppressed 



non-renormalizable terms, they can in general give rise to a large tadpole term in the 
potential proportional to the heavy scale: M 2 usy M P (S + S*) . Technically, the tadpole 
is generated through higher-order loop diagrams in which the non-renormalizable in- 
teractions participate as vertices together with the renormalizable ones. A discrete 
global symmetry like the one discussed above would forbid this term but would lead to 
the appearance of disastrous domain walls upon its unavoidable spontaneous break- 
down. The generated large tadpole reintroduces the hierarchy problem, since due to 
its presence the singlet v.e.v. gets a value (S) 2 ~ M susy M P . It appears that N = 1 
supergravity, spontaneously broken by a set of hidden sector fields, is the natural set- 
ting to study the generation of the destabilising tadpoles. A thorough analysis carried 
out in Ref. [JHj shows that the only harmful non-renormalizable interactions are ei- 
ther even superpotential terms or odd Kdhler potential ones. In addition, operators 
with more than six powers of the cut-off in the denominator are harmless. Finally, a 
tadpole diagram is divergent only if it contains an odd number of 'dangerous' vertices. 
The solution of the /z-problem in the framework of the NMSSM could be rendered 
a viable one if the destabilization problem were circumvented. What is needed is a 
suitable symmetry that forbids the dangerous non-renormalizable terms and allows 
only for tadpoles of order Mj ' usy (S + S*). This symmetry should at the same time 
allow for a large enough -Z^-breaking term in the scalar potential in order to destroy 



the unwanted domain walls |lT |. 



An alternative approach is to impose a symmetry which, although it does not 
forbid the dangerous non-renormalizable terms, only allows for higher-order tadpole 
graphs that give a n-loop-suppressed term (w ^ 2)n M 2 usy M P (S + S*). A case of partic- 
ular interest is when the cubic self-interaction for the singlet in eq. ([[]) is forbidden by 
the symmetry. Actually, it should be noted that if the underlying theory is a Grand 
Unified one (GUT), although a candidate for the singlet exists, a cubic term does 
not ariseQ- On the other hand, this case is truly minimal in the sense that, apart 
from promoting the /i-parameter into a field, no new renormalizable terms appear in 
the superpotential. Of course, a substitute is needed for the twofold role played by 
the cubic term, namely, its contribution to the mechanism generating the v.e.v. of S 
through the soft susy breaking terms and the breaking of the Peccei-Quinn symmetry 
present when k = 0. This role can be played by the tadpole. (Note that this is not 
included in the k — > limit of the existing NMSSM analyses, which up to now have 
ignored the tadpole term |5|, |6|, 12].) Recently, a viable solution along these lines was 



proposed based on discrete R-symmetries [[HJ . The renormalizable superpotential for 
this new minimal supersymmetric extension of the Standard Model (nMSSM) is given 
by 

W = XSH X H 2 + Y U QU C H X + Y d QD c H 2 + Y e LE c H 2 . (2) 

Apart from the usual Baryon and Lepton number, it possesses two additional global 
continuous symmetries, namely, an anomalous Peccei-Quinn symmetry U(1) p q with 

^dn Eq for example, matter and Higgs fields are contained in the 27 representation together with 
a singlet. Although the standard trilinear singlet-Higgses term is present in the 27 3 coupling, no 
singlet cubic term arises. The same is true for Eq embedings of 50(10) and SU(5). 



charges 

Q(-l), f/ c (0), D c (0), L(-l), £ c (0), ^(1), F 2 (l), 5(-2) (3) 

and a non- anomalous i?-symmetry U(1)r with charges 

Q(l), t/ c (l), D c (l), L(l), £ c (l), HtiO), H 2 (0), 5(2). (4) 

One of the solutions worked out consists in imposing the discrete sub-symmetry 2Z 5 p 
of the U(1)r' combination R' = 3R + PQ on the complete theory, including non- 
renormalizable operators. The charges under 2Z§r are 



(H h H 2 ) - 


-> a(ifi,# 2 ), 


(Q,L) - 


- « 2 (Q,L), 


(U C ,D C ,E C ) - 


-> a 3 (U c ,D c ,E 


S - 


-> a 4 ^, 


w - 


■+ aW, 



(5) 



where a = e 2l7r / 5 . An adequately suppressed linear term is generated at six- loop level 
by combining the non-renormalizable Kahler potential terms \\S 2 HiH 2 /Mp + h.c. 
and X 2 S (HiH 2 ) 3 / Mp + h.c. with the renormalizable superpotential term \SHxH 2 : 

V tadpole ~ -^—AiAaA 4 ^^,^^^ + 5*). (6) 

This tadpole has the desired order of magnitude 0(M susy ) if A1A2A 4 ~ 10 -3 . 

The goal of this paper is the phenomenological study of this nMSSM, where a 
gauge singlet superfield is added to the MSSM spectrum and a global Z 5R symmetry 
is imposed on the complete theory, resulting in the superpotential of eq. (0) and the 
tadpole term of eq. (|6]). In section 2 we review the general properties of the parameter 
space of the model. Phenomenological aspects of the nMSSM are addressed in section 
3. The neutralino sector, including the (quasi-pure) singlino, is studied in some 
detail. Also bounds on CP-even Higgs masses versus their couplings to gauge bosons 
are displayed along with Higgs production cross sections at future hadron colliders. 
Section 4 contains our main conclusions. 



2 Model set-up 

The tree-level Higgs scalar potential, namely, the potential which contains the scalar 
fields H x = (#?, Hi), H 2 = (# 2 + , # 2 °) and S, has the form : 



V® = V F + V D + V soft + V tadpole , (7) 

V F = \\\ 2 \(\Hi\ 2 + \H 2 \A\S\ 2 + \Hi\ 2 \H 2 \ 2 



\\\'lH»H»HrH+ + h.c) (8) 



,2 _l „2 



T/ Pi + ff 2 
V^D = o 



l n l 2 l w l 2 
-Hi — \rl 2 \ 



12 ^2 



ffl^^l 2 , (9) 



Vsoft = m 2 Hl \H l \ 2 + mj l2 \H 2 \ 2 + m 2 s \S\ 2 

+ f XA. X SH X H 2 + h.c). (10) 



In what follows, we shall assume a phenomenological point of view and write the 
generated tadpole as 

v tadpole = e(s+s*), (ii) 

where £ is treated as a free parameter. 

In order to obtain the correct upper limits on the Higgs boson masses (cf. section 
|3.2| ) radiative corrections to the tree-level potential have to be considered. Let us 
introduce a scale Q ~ M susy and assume that quantum corrections involving momenta 
P 2 ~ Q 2 have been evaluated, e.g. by the integration of the Renormalization Group 
Equations (RGEs) of the parameters from initial values at the GUT scale down to 
the scale Q. One is then left with the computation of quantum corrections involving 
momenta p 2 < Q 2 . The effective potential V e ff can be developed in powers of % or 
loops as 

V eff = V<® + V^ + V^ + .... (12) 

The tree-level potential V^ is given by eq. ([?]). The one-loop corrections to the 
effective potential read as 



V {1) = STrM 4 

64tt 2 




(13) 



where M 2 is the field dependent squared mass matrix (in our analysis, we take only 
top/stop loops into account). Next, we consider the dominant two-loop corrections. 
These will be numerically important only for large susy breaking terms compared to 
the Higgs v.e.v.'s hi, hence we can expand in powers of hi. Since the terms quadratic 
in hi can be absorbed into the tree-level soft terms, we just consider the quartic 
terms, and here only those which are proportional to large couplings: terms ~ a s hf 
and ~ h\. Finally, taking only leading logarithms (LLs) into account, the expression 
for V^ reads 

V% = z(^h) 2 hi(v™ S -~hi)t 2 , (14) 



x 16tt 2 / z V 2 

where t = In ( ^- J , rrit being the top quark mass. One- loop corrections to the tree- level 
relations between bare parameters and physical observables, once reinserted in the 
one-loop effective potential, also appear as two-loop effects. These are: corrections to 
the kinetic terms of the Higgs bosons, which lead to a wave function renormalization 
factor Zh 2 m front of the D^HiD^Hi term with, to order h 2 



^ 2 = 1 + 3 16^ t; (15) 

and corrections to the top quark Yukawa coupling with, to orders h\, a s 

h t (m t ) = h t (Q) (1 + 3-^ (32vra s - % - hfj t) . (16) 

In general, the parameters A, A\ and £ could be complex. However, by redefining 
the fields H 2 (or Hi) and S, one can always get - without loss of generality - that 
\A\,^ 3 G R. Note that in the NMSSM with the cubic singlet superpotential term 
^kS 3 one has to further assume that the combination \k* (or, equivalently, A\/A K ) 
is real || ||. By SU(2)l X U(l)y gauge invariance one can get rid of the phase of Hi, 
by taking {H{ ) = and hi = (H® ) G R + . One can then show that the condition for a 
local minimum with (H 2 ) = is equivalent to a positive mass squared for the charged 
Higgs. It has been proven that a sufficient condition is A < g 2 [|nj which, as we shall 



see below, is always verified in the universal case. By taking h 2 = (#2) = P2^ 2 , 
s = (S) = poe 1 ^ and minimizing the complete (two-loop) effective potential with 
respect to <fto arid <p 2 , we find that there is one and only one global vacuum for 
which the two phases relax to zero, i.e., 0o = <ft 2 = 0. This implies that there is no 
spontaneous CP -violation. Therefore one can choose hi G R + and h 2 ,s G R. This 
result distinguishes the nMSSM from the usual NMSSM where loop corrections can 



generate spontaneous CP-violation |fi5 



The soft terms of the model can be constrained by requiring universality at the 
GUT scale. The independent parameters of the model are then a universal gaugino 
mass Mi/2 (always positive in our convention), a universal mass for the scalars m\, 
a universal trilinear coupling A (either positive or negative), the (positive) Yukawa 
coupling Ao at the scale Mqut and the tadpole coefficient £. The (well-known) value of 
the Z-boson mass fixes one of these parameters with respect to the others, so that we 
end up with four free parameters at the GUT scale, i.e., as many as in the MSSM with 
universal soft terms. In principle, one could choose the same set of free parameters as 
in the MSSM, i.e., Mi/ 2 , m^, A and tan/3(= j^-), with A, s and £ being determined 
by the three minimization equations, demanding also radiative electroweak symmetry 
breaking 0. However, this appears to be a non-trivial issue, as A also influences the 
running of the RGEs of the soft parameters between the GUT and the EW scale. In 
other terms, one would need a lot of fine tuning of the dimensionful A Q in order to 
get a dimensionless parameter, A, of the desired value at the EW scale. Therefore, in 
the case of universality, we conveniently adopt in our numerical analysis the following 
input parameters: ml/Mi/ 2 , A /Mi/ 2 , i/My 2 and A (tan/5 and s being calculated 
from the minimization of the potential and the overall scale Mi/ 2 fixed by Mz)- 

If one requires the absence of a Landau singularity for A below the GUT scale, 
one obtains an upper bound on A at the EW scale. This upper bound depends on 
the value of the top quark Yukawa coupling h t , i.e., on tan/3 (cf. fig. 1). Requiring 
furthermore, universality at the GUT scale, one ends up with the more restrictive 
constraint A < .3, higher values leading to unphysical global minima of the effective 
potential. 
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Figure 1: Upper bound on A as a function of tan/? for m? t ° e 
The width of the curve is due to the uncertainty on mf° e . 



173.8 ± 5.2 GeV JL8 



Let us now briefly address the problem of Charge and Color Breaking (CCB) min- 
ima. The most dangerous CCB direction involves the trilinear coupling h e A e E^iLiHi 
where h e denotes the electron Yukawa coupling (~ 10~ 5 ), E R i is the right-handed 
selectron, L x the left-handed slepton doublet of the first generation, and Hi the cor- 
responding Higgs doublet. From the absence of a non-trivial minimum of the scalar 
potential in the .D-flat direction \Er^\ = \L\\ = |i?2|, the following inequality among 
the soft susy breaking terms can be derived [[§]: 



Ai<3 



m F + m T + m 



;i7) 



where m%, m\ and m\ are the soft susy breaking mass terms associated with the three 
fields above. If the inequality ( |TTD is violated, the fields develop v.e.v.'s of O (A e /h e ) 
and the depth of the minimum is of O (A^/h 2 e ). Accordingly, ([H]) has to be imposed 
at a scale Q ~ A e /h e ~ 10 7 GeV. Assuming universal soft terms at the GUT scale, 
(|T7l) then becomes 0: 



(A) - 0.5M 1/2 ) 2 < 9m 2 + 2.67 M{ 



/2- \ LO ) 

So-called UFB directions (for Unbounded From Below, which actually never occurs in 
the universal case) can also be considered. Assuming universality for the soft terms, 



the absence of a global minimum in these directions typically implies [|16" , |17|1 



m 

M 1/2 



>1. 



(19) 



However, the tunnelling rate from the standard EW minimum to a UFB one is in 
general quite small fl6f , so that this constraint can be avoided if one is ready to 



assume that the standard EW vacuum is metastable. 



3 Phenomenological aspects of the nMSSM 

3.1 Singlino LSP and additional cascades 

The nMSSM contains additional gauge singlet states in the Higgs sector (one neutral 
CP-even and one CP-odd state) and in the neutralino sector (a two component Weyl 
fermion). These states are mixed with the corresponding ones of the MSSM, and the 
physical states have to be obtained from the diagonalization of the mass matrices in 
each sector. In the basis (B, W 3 , H®, H®, S) the (symmetric) neutralino mass matrix 
reads as 



/ Mi 



M ( 



\ 



-gihi/y/2 


gih 2 /V2 


W 2 92hi/y/2 


-g 2 h 2 /V2 





As Xh 2 




Xhi 








\ 



/ 



(20) 



Note that the diagonal singlino (S) mass term is zero. Furthermore, the singlino 
mixings (with Higgsinos) are proportional to A, which, as remarked earlier, turns out 
to be quite small, especially in the universality scenario. Consequently, the singlino 
state appears to be an almost pure singlet state with a very small mass, so that it 
is always the lightest supersymmetric particle (LSP) of the model. Actually, in the 



universal case, we find: few MeV < 



< 



rv ^g rv< 



3 GeV, with a singlet component 



> 



o. 



This state has only small couplings to the gauge bosons and to the other sparticles. 
(We have explicitly checked that its contribution to the invisible Z-boson width is 
< 4.2 MeV [jj~8|.) Thus, the production cross sections of the singlino are small and it 
seems to be nearly impossible to observe this particle in any experiment, this rendering 
the nMSSM apparently similar to the ordinary MSSM. However, as the singlino is the 
LSP of the model, it will appear at the end of all sparticle decay chains, giving rise 
to additional cascades compared with the MSSM signals. Such additional cascades 
are common to many supersymmetric extensions involving singlets |i~2 , |l9fl . It should 
be noticed, however, that unlike in the NMSSM, where the singlino LSP scenario 
requires strong constraints on the parameter space (i.e., M\/ 2 ^> m ,A ) fl2|], the 
singlino is always the LSP in the nMSSM. 

Besides, by assuming universality at the GUT scale, we find that the next-to- 
lightest supersymmetric particle (NLSP) is always the second lightest neutralino, 
which turns out to be a quasi-pure bino (B). Depending on the region of the param- 



eter space under scrutiny, the following channels can play a role in the NLSP — > LSP 
cascade: 

• B —* Svv (sneutrino/Z exchange) giving an invisible cascade; 

• B —* Sl + l~ (slepton/Z exchange) where the leptons could be mainly r's, the 
stau being lighter than the other sleptons; 

• B — > Sqq (squark/Z exchange) the branching ratio being quite small ( < 10%), 
as the squarks are usually heavy; 

• B — > SZ if the B is heavy enough; 

• B — > SS where S is a light quasi-pure singlet Higgs boson, decaying to bb or tt 
depending on its mass; 

• B — > S7 through loops. 



The properties of these cascades have been analyzed in detail in Ref. |12| for the case 
of the NMSSM and most of the results can equally apply to the case of the nM SSM. 
As for experimental searches, high multiplicity events have been under investigation 
already in the context of models with gauge mediated supersymmetry breaking P(J 
or with i?-parity violation ||21|| . In principle, small A's ( <• 1CT 4 ) could give rise to a 



delayed NLSP — > LSP transition, i.e. a displaced neutral vertex [|T2| , [19|j . However, 
such values of A are disfavoured if one wants the tadpole term of eq. (|6]) to be large 
enough, so that displaced neutral vertices are not expected as typical signatures of 
the nMSSM. 

3.2 Higgs couplings and mass bounds 

The Higgs sector of the nMSSM consists of three CP-even neutral states, denoted 
by Si with masses ra\ < m 2 < m 3 , plus two CP-odd neutral states, labelled as Pi 
with masses m! x < m' 2 . The tree-level mass matrix for the CP-even states in the basis 
(ReH°, ReH$, ReS) reads 



g^j-XsAxtanP {2\ 2 - g 2 )h 1 h 2 + XsA x A(2As/i! + A x h 2 ) \ 
s = I g 2 h\ - XsA x cot j3 X(2Xsh 2 + A x h ± 



-^^ - a£ J 



(21) 



where g 2 = (g 2 + g 2 )/2. In the reminder of this section, we study the upper bounds on 
the lightest CP-even states with general soft susy breaking terms, in a non-universal 
scenario. By taking into account the full one-loop and the dominant two-loop top/stop 
corrections displayed in section 2, and assuming hi <C M susy , one obtains the following 
upper limit on the lightest CP-even Higgs mass: 



ml < M\\ cos 2 2/3 + 



2A 2 



where 



gl + 92 
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Figure 2: Upper bound on m\ [GeV] versus tan (5 for mf° e = 173.8±5.2 GeV (straight 
dashed, dotted line respectively), and M susy < 1 TeV. 



The only difference between the MSSM bound [[22] and eq. (j22|) is the 'tree-level' 
contribution ~ A 2 sin 2 2/5. This term is important for moderate values of tan (3. Hence, 
the maximum of the lightest Higgs mass in the nMSSM is not obtained for large 
tan/3 values, like in the MSSM, rather for moderate ones (cf. fig. 2). In contrast, the 
radiative corrections are identical in the nMSSM and in the MSSM. In particular, 
the linear dependence in X t is the same in both models. Hence, from eq. (^), the 
upper bound on m 2 is maximised for X t = 6 (corresponding to A t = y/6M susy , the 



'maximal mixing' case) and minimized for X t — (corresponding to A t — 0, the 'no 
mixing' case). 

However, the upper limit on mi is not necessarily physically relevant, since the 
coupling of the lightest CP-even Higgs boson to the Z-boson can be very small. 
Actually, this phenomenon can also appear in the MSSM, if sin 2 (/3 — a) is small. In 
this case though, the CP-odd Higgs boson A is necessarily light {tua ~ m/, < Mz at 
tree-level) and the process e + e" — ► Z — ► /iA can be used to cover this region of the 
MSSM parameter space. In the nMSSM, a small gauge boson coupling of the lightest 
Higgs S\ is usually related to a large singlet component, in which case no (strongly 
coupled) light CP-odd Higgs boson is available. Hence, Higgs searches in the nMSSM 
have to possibly rely on the search for the second lightest Higgs scalar S 2 . 

Let us define the reduced coupling Ri as the square of the ZZSi coupling divided 
by the corresponding standard model Higgs - Z boson coupling 

Ri = (S a cos/? + S l2 sin/?) 2 , (25) 

where Sn, Si 2 ar e the Hi, H 2 components of the CP-even Higgs boson Si, respectively. 
Evidently, we have < Ri < 1 and unitarity implies 

£fl* = 1 - (26) 

We are interested in upper bounds on the two lightest CP-even Higgs bosons Si^ 2 . 
These are obtained in the limit where the third Higgs, S3, is heavy and decoupled, 
i.e., i?3 ~ (this scenario is similar the so-called decoupling limit in the MSSM: the 
upper bound on the lightest Higgs h is saturated when the second Higgs H is heavy 
and decouples from the gauge bosons). In this limit, we have Ri + R 2 — 1. 

In the regime Ri > 1/2, experiments will evidently first discover the lightest Higgs 
(with mi < 133.5 GeV for mf oe = 173.8 GeV and M susy = 1 TeV). The 'worst case 
scenario' in this regime corresponds to mi ~ 133.5 GeV and Ri ~ 1/2: the presence 
of a Higgs boson with these properties has to be excluded in order to test this part 
of the parameter space of the nMSSM. 

In the regime Ri < 1/2 (i.e. 1/2 < R 2 < 1) the lightest Higgs may escape 
detection because of its small coupling to gauge bosons, and it may be easier to look 
for the second lightest Higgs S 2 . In fig. 3 we show the upper limit on m 2 as a function 
of R 2 as a thin straight line. For R 2 — > 1 (i.e. Ri — > 0), the upper limit on m 2 
is actually given by the previous upper limit on mi, even if the corresponding Higgs 
boson is the second lightest one. For R 2 — > 1/2, on the other hand, m 2 can be as large 
as 190 GeV. However, one finds that the upper limit on m 2 is saturated when the 
mass mi of the lightest Higgs boson tends to 0. Clearly, one has to take into account 
the constraints from Higgs boson searches which apply to reduced couplings R < 1/2, 
i.e., lower limits on mi as a function of Ri ~ 1 — R 2 , in order to obtain realistic 
upper limits on m 2 versus R 2 . The dotted curves in fig. 3 show the upper limit on 
m 2 as a function of R 2 for different fixed values of mi (as indicated on each curve). 
They can be used to obtain upper limits on the mass m 2 , in the regime Ri < 1/2, for 
arbitrary experimental lower limits on the mass m 1 versus R-\_. For each value of the 

10 



coupling Ri, which would correspond to a vertical line in fig. 3, one has to find the 
point where this vertical line crosses the dotted curve associated to the corresponding 
experimental lower limit on mi. Joining these points by a curve leads to the upper 
limit on 7772 as a function of R2. We have indicated as a thick straight line in fig. 3 



the present LEP-II bounds pUfl , which give, in the 'worst case' scenario, an upper 
limit on m 2 of ~ 160 GeV for R 2 — 1/2. 




Figure 3: Upper limits on the mass mi (in GeV) against R2, for different values of 
mi (as indicated on each line in GeV), assuming vnP t oe = 173.8 GeV and M susy < 1 
TeV. R\ — 1 — R2 is shown on the the top axis. The thick straight line corresponds 
to LEP-II lower limits on mi vs. R\. 



Lower experimental limits on a Higgs boson with R > 1/2 restrict the allowed 
regime for m 2 (for R 2 > 1/2) in fig. 3 from below. The present lower limits on m 2 
from LEP-II are not visible in fig. 3, since we have only shown the range m 2 > 130 
GeV. Possibly Higgs searches at Tevatron Run II will push the lower limits on 777-2 
upwards into this range. This would be necessary if one aims at an exclusion of 
this regime of the nMSSM. Then, lower limits on the mass 7712 - for any value of R2 
between 1/2 and 1 - of at least 133.5 GeV are required. The precise experimental 
lower limits on rri2 as a function of R2, which would be needed to this end, will depend 
on the achieved lower limits on m,\ as a function of R\ in the regime R\ < 1/2. 

In principle, from eq. (|26"1), one could have R2 > Ri with R 2 as small as 1/3. 
However, in the regime 1/3 < R2 < 1/2, the upper bound on 7772 as a function of 
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i?2 for different fixed values of m>i can only be saturated if R\ = R2. It is then 
sufficient to look for the lightest Higgs S\ (i.e. for a Higgs boson with a coupling 
1/3 < R < 1/2 and a mass m <• 133.5 GeV) to cover this region of the parameter 
space of the nMSSM. 

One can notice that these results are the same as those obtained in the NMSSM 
p4| . This comes from the fact that the non-singlet part of the CP-even mass matrix 
(^T|) as well as the singlet /non-singlet mixing terms are the same in both models, 
giving the same upper limit on mi, eq. ( ^2]) and fig. 2. Even though they are not 
similar, the CP-even singlet mass term (the 3x3 element in eq. fl2lD ) is 'free' in 
both models, i.e., it can take any value between and ~ 1 TeV. This explains why 
the curves displayed in fig. 3 are the same in both models, the upper bounds on m 2 
steaming from degenerate singlet /non- singlet mass terms, degeneration lifted by the 
mixings (off-diagonal 1x3 and 2x3 terms in eq. (p2|)). 

Thus, the phenomenological potential of the nMSSM is not less exciting in the 
Higgs sector than it is in the sparticle sector. On the one hand, the necessary (but 
not sufficient) condition for testing the complete parameter space of the nMSSM is 
to rule out a CP-even Higgs boson with a coupling 1/3 < R < 1 and a mass below 
135 GeV. On the other hand, the sufficient condition (i.e., the precise upper bound 
on 1712 versus R2) depends on the achieved lower bound on the mass of a 'weakly' 
coupled Higgs (with < R < 1/2) and can be obtained from fig. 3. At the Tevatron 
Run II this would probably require an integrated luminosity of up to 30 fb _1 p5] . If 
this cannot be achieved, one has to wait for the advent of the LHC, in order to know 
whether the nMSSM is actually realized in nature. 

Imposing universality of the soft terms at the GUT scale, CCB constraint (eq. fll7|)) 
and LEPII experimental constraints on Higgs masses |23[] one finds the more restric- 
tive bound mi < 122 GeV under the same hypotheses (mf oe = 173.8 ± 5.2 GeV and 
M susy < 1 TeV) and m 2 < 135 GeV in the case where Si is mainly singlet (i?i < R2). 

An intriguing example of possible evidence of the nMSSM at future hadron-hadron 
colliders is the following. Recall that a crucial signature for an intermediate Higgs bo- 
son at both the Tevatron and the LHC is the one produced via qq' — > W^ Higgs (with 
a smaller contribution from qq — ► Z Higgs as well), with the gauge vector yielding 
high transverse-momentum and isolated leptons and the Higgs scalar decaying into bb 
pairs. Now, let us imagine that Higgs searches in this mode have finally revealed the 
evidence of a light scalar Higgs resonance, but no further Higgs states are detected 
up to well above the EW scale. (In this respect, the LHC is a better example to 
illustrate, as compared to the Tevatron, because of its much extended scope in mass.) 
This scenario could be realized in the MSSM, if the latter is in the above mentioned 
decoupling regime, where one has vnn± « m# ~ rriA ^> Mz and the strength of 
the lightest Higgs boson couplings to the gauge vector bosons W ± and Z approaches 
unity, Rh = sin 2 (/5 — a) ~ 1. Indeed, this phenomenology can also be realized in the 
nMSSM when, e.g., R\ — 1 and R 2 = OR. Under these circumstances, the mass of the 



2 Further notice that, no matter their actual mass value, pseudoscalar neutral Higgs states of 
either model cannot be produced via the above two processes, as their coupling to W^ and Z 
vectors is prohibited at tree-level. 
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h scalar would only depend on tan /3 and the minimum of the Higgs production cross 
section is obtained, in both susy models, in correspondence of the maximum Higgs 
mass. Fig. 4 shows this dependence for both the MSSM and the nMSSM at the LHC, 
with \/s = 14 TeV. (The trend is very similar at the Tevatron, where y/s = 2 TeV.) 
Over a broad range in tan /3, the production rates of the latter are significantly below 
those of the former, with a minimum at tan/3 ~ 2.7 (for m t = 173.8 GeV), owning to 
the peculiar tan/3 dependence of the maximum value of m l5 as seen in fig. 2 (recall 
instead that the maximum mass of the lightest Higgs boson of the MSSM increases 
monotonically with tan (3 |22| , hence the cross section here decreases correspondingly, 
for tan (3 < 40). 



P P 



V Si (V=W,Z) at NLO in the newMSSM 
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Figure 4: Minimum of the production cross section of a neutral scalar Higgs with 
SM-like couplings to W^ and Z gauge bosons at the LHC, as a function of tan/3, 
in correspondence of the maximum values of its mass, in both the MSSM and the 
nMSSM. 

Thus, if tan/? has already been measured in another context (e.g., in susy sparti- 
cle processes), to detect an isolated scalar resonance in the bb channel (the dominant 
decay model in either model for Mmggs < 133.5 GeV), at a rate well below the mini- 
mum one predicted by the MSSM, could imply that physics beyond the latter would 
be realized in nature. In fact, the reader should recall that the BR(Higgs — > bb) is 
basically the same in both models |26| and the combined uncertainties on the rates 
of the latter and of the production modes (due to higher-order effects, parton distri- 
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bution functions, hard scale dependence, etc.) are of the order of just a few percent 



27 1 . In contrast, the differences between the MSSM and nMSSM rates can be as 



large as a factor of two, in the vicinity of tan/5 ~ 2.7, and well above the mentioned 
uncertainties for tan/3 up to 10 or so. Besides, kinematical analysis of the bb system 
might provide further evidence in this respect, if the mass resolution is larger than 
the difference between the Higgs mass values, as predicted by the two models for a 
given tan f3. Similar arguments can be made for the case of the gg — ► Higgs — ► 77 
signature too, however, the production and decay phenomenology is here much more 
involved (because loop processes take place at either stage), so that we leave it aside 
for future consideration [^ . 



4 Summary and conclusions 

We have studied the phenomenology of the nMSSM, which promotes the /i-parame- 
ter into a singlet superfield, hence solving the well-known /i-problem of the MSSM. 
Besides, cubic self-couplings and possible dimensionful couplings are avoided in the 
nMSSM thanks to a global discrete R-symmetry, in turn broken by supergravity- 
induced tadpole corrections, which solves both the so-called 'domain wall' and 'axion' 
problems. The new model is truly minimal, in the sense that - despite incorporating 
new fields - it can be parametrised by the same number of inputs as in the MSSM in 
the universal case. Its phenomenology stands out quite different from that of both the 
ordinary MSSM and the so-called Next-to-Minimal Supersymmetric Standard Model 
(NMSSM) - where the cubic self-interacting coupling is instead present. In particular, 
the following aspects emerged as crucial from our analysis. 

• Assuming an accuracy up to the dominant top-stop contributions at two-loop 
level and depending on the (reduced) couplings of the two lighter CP-even Higgs 
bosons to the Z-boson, we have found that the upper limit on the mass of the 
lightest Higgs state of the nMSSM can be 133.5 GeV, in correspondence of the 
central value of the top mass (i.e., m t = 173.8 GeV), that is, about 10 GeV 
higher than the MSSM value and within the Tevatron Run-II reach. In addition, 
the upper limit on the mass of the lightest 'observable' Higgs boson (i.e., the 
next-to-lightest one, when the lightest one couples invisibly to the Z-boson) 
could be as high as 160 GeV but still within the LHC scope. 

• To remain with the Higgs sector, we also have described a benchmark example 
that could allow one to phenomenologically distinguish the nMSSM from the 
MSSM in the search for the lightest Higgs state at future hadron colliders, 
such as the Tevatron Run-II and the LHC. If only one neutral Higgs state is 
accessible through associated production with an EW gauge vector W ± or Z, 
via its bb decays, the knowledge of tan/3 and of the production rate of such 
Higgs process could be enough to assign such a Higgs state to one or the other 
of the two models, even prior to the investigation of the mass resonance that can 
eventually be reconstructed from the bb system. Production and decay studies 
of all other Higgs states of the nMSSM are now in progress [2j| • 
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• Despite the remarkable dissimilarities seen so far between the nMSSM and the 
MSSM, one might quite rightly question that the phenomenology of the nMSSM 
is very similar to that of the NMSSM, as far as the Higgs sector is concerned. 
However, a dramatic difference is revealed between these two models, if one 
investigates CP-violation effects. In fact, no matter the actual value of A, a 
peculiar feature of the nMSSM is the following: Contrary to the NMSSM case, 
spontaneous CP-violation cannot occur in the Higgs sector of the nMSSM, nei- 
ther at tree-level nor at one-loop. 

• Finally, further differences between the nMSSM and any other model can be 
appreciated in the sparticle sector. In fact, here, two concurrent aspects ren- 
der the phenomenology of the former both more 'spectacular' and 'natural', 
in comparison to the MSSM and the NMSSM, respectively. Firstly, the light- 
est neutralino appears to be an almost pure singlet state. Secondly, such a 
state (the 'singlino') is always the LSP of the theory, with a very small mass 
(varying from a few MeV to a few GeV). The consequence is twofold. On the 
one hand, in comparison to the NMSSM, the singlino LSP scenario requires no 
strong constraints to be imposed on the parameter space. On the other hand, 
in comparison to the MSSM, any sparticle decay chain involves a further step, 
the NLSP — > LSP transition, giving rise to additional cascades. 

• The stimulating issue that such an LSP could be a good dark matter candidate 
or, alternatively, could be excluded from cosmological arguments, also deserves 
attention. However, a quantitative analysis in this respect was far beyond the 
intention of this note. 

Concluding, the nMSSM is, at the same time, the best theoretically motivated and the 
most economical susy extension of the SM. Here, we have pointed out differences or 
similarities between the new model, the NMSSM and the MSSM, in both the Higgs 
and neutralino sectors. Future collider experiments, at Tevatron Run-II and LHC, 
will be able to prove whether or not the nMSSM is the realization of Supersymmetry 
that nature has chosen. 

Note While finalising this paper, whose main results were made public already in 



Ref. p9|| , we became aware of Ref. [pO] , by C. Panagiotakopoulos and A. Pilaftsis, 
which deals with a similar subject. In this respect, although we agree with the more 
general statements given in this other paper, we would like to remark that there only 
the top-stop one-loop corrections were used, whereas here we have calculated also 
the corresponding two-loop contributions. This explains the significant discrepancy 
between their and our upper limit on the lightest Higgs boson mass: mi < 150 GeV 
versus mi < 133.5 GeV. 
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